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ABSTRACT 

Mergers of black-hole binaries are expected to release large amounts of energy 
in the form of gravitational radiation. However, binary evolution models predict 
merger rates too low to be of observational interest. In this paper we explore 
the possibility that black holes become members of close binaries via dynamical 
interactions with other stars in dense stellar systems. In star clusters, black 
holes become the most massive objects within a few tens of millions of years; 
dynamical relaxation then causes them to sink to the cluster core, where 
they form binaries. These black-hole binaries become more tightly bound 
by superelastic encounters with other cluster members, and are ultimately 
ejected from the cluster. The majority of escaping black-hole binaries have 
orbital periods short enough and eccentricities high enough that the emission 
of gravitational radiation causes them to coalesce within a few billion years. 
We predict a black-hole merger rate of about 1.6 x 10~ 7 per year per cubic 
megaparsec, implying gravity wave detection rates substantially greater than the 
corresponding rates from neutron star mergers. For the first generation Laser 
Interferometer Gravitational- Wave Observatory (LIGO-I), we expect about one 
detection during the first two years of operation. For its successor LIGO-II, the 
rate rises to roughly one detection per day. The uncertainties in these numbers 
are large. Event rates may drop by about an order of magnitude if the most 
massive clusters eject their black hole binaries early in their evolution. 



1. Introduction 

The search for gravitational waves will begin in earnest in January 2002, when LIGO-I 
(Abramovici et al. 1992) becomes fully operational (K. Thorne, private communication). 
The appearance of this new and wholly unexplored observational window challenges 
physicists and astronomers to predict detection rates and source characteristics. Mergers of 
neutron-star binaries are widely regarded as the most promising sources of gravitational 
radiation, and estimates of neutron star merger rates (per unit volume) range from 
11 ~ 1.9 X 10" 7 h 3 yr" 1 Mpc~ 3 (Narayan et al. 1991; Phinney 1991; Portegies Zwart 
& Spreeuw 1996) , where h = iJ /100kms _1 Mpc -1 , to roughly ten times this value 
(Tutukov & Yungelson 1993; Lipunov et al., 1997). However, even with the most optimistic 
assumptions, we can expect a LIGO-I detection rate of only a few neutron star events per 
millennium. 

Inspiral and merger of black-hole binaries are considerably more energetic events than 
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neutron star mergers, due to the higher masses of the black holes (Tutukov & Yungelson 
1993; Lipunov et al. 1997). Black-hole binaries can result from the evolution of two 
stars which are born in a close binary, experience several phases of mass transfer, and 
subsequently survive two supernovae (Tutukov & Yungelson 1993). Calculations of event 
rates from such field binaries depend sensitively on many unknown parameters and much 
poorly understood physics, but the models generally predict a black hole merger rate 
1Z rS 2 x 10~ 9 h 3 yr^ 1 Mpc~ (Tutukov & Yungelson 1993; Portegies Zwart & Yungelson 
1998; Bethe & Brown 1999), substantially lower than the rate for neutron stars. An 
alternative possibility, which we explore here, is that black holes become members of close 
binaries not through internal binary evolution, but rather via dynamical interactions with 
other stars in a dense stellar system. 



2. Black hole binaries in star clusters 

Black holes are the products of stars with initial masses exceeding ~2O-25M (Maeder 
1992; Portegies Zwart et al. 1997). A Scalo (1986) mass distribution with a lower mass limit 
of 0.1 M and an upper limit of 100 M Q has 0.071% of its stars more massive than 20 M , 
and 0.045% more massive than 25 M . A star cluster containing N stars thus produces 
~ 6 x 10~ 4 iV black holes. Known Galactic black holes have masses m^ between 6M and 
18 M (Cowlay 1992). For definiteness, we adopt rabh = 10 M . 



2.1. Binary formation and dynamical evolution 

A black hole is formed in a supernova explosion. If the progenitor is a single star (i.e. 
not a member of a binary), the black hole experiences little or no recoil and remains a 
member of the parent cluster (White & van Paradijs 1996). If the progenitor is a member of 
a binary, mass loss during the supernova may eject the binary from the cluster potential via 
the Blaauw mechanism (Blaauw 1962), where conservation of momentum causes recoil in a 
binary which loses mass impulsively from one component. We estimate that no more than 
~ 10% of black holes are ejected from the cluster immediately following their formation. 

After ~ 40 Myr the last supernova has occurred, the mean mass of the cluster stars is 
(m) ~ 0.56 M (Scalo 1986), and black holes are by far the most massive objects in the 
system. Mass segregation causes the black holes to sink to the cluster core in a fraction 
~ (m)/mbh of the half-mass relaxation time. For a typical globular cluster, the relaxation 
time is ~ 1 Gyr; for a young populous cluster, such as R 136 (NGC 2070) in the 30 Doradus 
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region of the Large Magellanic Cloud (Massey & Hunter 1998), it is ~ 10 Myr. 

By the time of the last supernova, stellar mass loss has also significantly diminished 
and the cluster core starts to contract, enhancing the formation of binaries by three-body 
interactions. Single black holes form binaries preferentially with other black holes (Kulkarni 
et al. 1992), while black holes born in binaries with a lower-mass stellar companion rapidly 
exchange the companion for another black hole. The result in all cases is a growing 
black-hole binary population in the cluster core. Once formed, the black-hole binaries 
become more tightly bound through superelastic encounters with other cluster members 
(Heggie 1975; Kulkarni et al. 1992; Sigurdsson & Hernquist 1993). On average, following 
each close binary-single black hole encounter, the binding energy of the binary increases 
by about 20% (Hut et al. 1992); roughly one third of this energy goes into binary recoil, 
assuming equal mass stars. The minimum binding energy of an escaping black-hole binary 
may then be estimated as 

E b , miQ ~36W ^kT, (1) 

(m) 

where |/cT is the mean stellar kinetic energy and Wq = (m)\(f)o\/kT is the dimensionless 
central potential of the cluster (King 1966). By the time the black holes are ejected, 
(m) ~ 0.4 M Q . Taking Wo ~ 5-10 as a representative range, we find -E^min ~ 5000-10000 
kT. 

We have tested and refined the above estimates by performing a series of iV-body 
simulations within the "Starlab" software environment (Portegies Zwart et al. 1999: see 
http: : /www. sns. ias.edu/~starlab), using the special-purpose computer GRAPE-4 to 
speed up the calculations (Makino et al. 1997). For most (seven) of our calculations we 
used 2048 equal-mass stars with 1% of them ten times more massive than the average; two 
calculations were performed with 4096 stars. One of the 4096-particle runs contained 0.5% 
black holes; the smaller black-hole fraction did not result in significantly different behavior. 
We also tested alternative initial configurations, starting some models with the black holes 
in primordial binaries with other black holes, or in primordial binaries with lower-mass 
stars. 

The results of our simulations may be summarized as follows. Of a total of 204 black 
holes, 62 (~ 30%) were ejected from the model clusters in the form of black-hole binaries. 
A total of 124 (~ 61%) black holes were ejected single, and one escaping black hole had 
a low-mass star as a companion. The remaining 17 (~ 8%) black holes were retained by 
their parent clusters. The binding energies Ef, of the ejected black-hole binaries ranged 
from about 1000 kT to 10000 kT in a distribution more or less flat in log Ej,, consistent with 
the assumptions made by Hut et al. (1992). The eccentricities e followed a roughly thermal 
distribution [p(e) ~ 2e], with high eccentricities slightly overrepresented. About half of the 
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black holes were ejected while the parent cluster still retained more than 90% ( £ 2 initial 
relaxation times) of its birth mass, and £ 90% of the black holes were ejected before the 
cluster had lost 30% (between 4 and 10 relaxation times) of its initial mass. These findings 
are in good agreement with previous estimates that black-hole binaries are ejected within 
a few Gyr, well before core collapse occurs (Kulkarni et al. 1993; Sigurdsson & Hernquist 
1993). 

We performed additional calculations incorporating a realistic (Scalo) mass function, 
the effects of stellar evolution, and the gravitational influence of the Galaxy. Our model 
clusters generally dissolved rather quickly (within a few hundred Myr) in the Galactic tidal 
field. Clusters which dissolved within ~ 40 Myr (before the last supernova) had no time to 
eject their black holes. However, those that survived beyond this time were generally able 
to eject at least one close black-hole binary before dissolution. 

Based on these considerations, we conservatively estimate the number of ejected 
black-hole binaries to be about 10~ 4 iV per star cluster, more or less independent of the 
cluster lifetime. 



2.2. Characteristics of the binary population 

The energy of an ejected binary and its orbital separation are coupled to the 
dynamical characteristics of the star cluster. For a cluster in virial equilibrium, we have 
kT = 2E'kin/3iV = —E pot /3N = GM 2 /6A^r vir , where M is the total cluster mass and r vir is 
the virial radius. A black-hole binary with semi-major axis a has Ef, = Gm1 h /2a, so 

^ = 3A^V^, (2) 

kT V M ) a ' V ; 

establishing the connection between a and the bulk parameters of the cluster. 

In computing the properties of the black-hole binaries resulting from cluster evolution in 
§0, it is convenient to distinguish three broad categories of dense stellar systems: (1) young 
populous clusters, (2) globular clusters, and (3) galactic nuclei. Table [I] lists characteristic 
parameters for each. The masses and virial radii of globular clusters are assumed to be 
distributed as independent Gaussians with means and dispersions as presented in the 
table; this assumption is supported by correlation studies (Djorgovski & Meylan 1994). 
Table [I] also presents estimates of the parameters of globular clusters at birth (bottom 
row), based on a recent parameter-space survey of cluster initial conditions (Takahashi & 
Portegies Zwart 2000); globular clusters which have survived for a Hubble time have lost 
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60% of their initial mass and have expanded by about a factor of three. We draw no 
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distinction between core-collapsed globular clusters (about 20% of the current population) 
and non-collapsed globulars — the present dynamical state of a cluster has little bearing on 
how black- hole binaries were formed and ejected during the first few Gyr of the cluster's 
life. 



3. Production of gravitational radiation 

An approximate formula for the merger time of two stars due to the emission of 
gravitational waves is given by Peters & Mathews 1963): 
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W g « 150 Myr — £ — (1 - e 2 ^ 2 . (3) 



The sixth column of Table |I| lists the fraction of black- hole binaries which merge within a 
Hubble time due to gravitational radiation, assuming that the binary binding energies are 
distributed flat in logE^ between 1000 kT and 10000 kT, that the eccentricities are thermal, 
independent of E^, and that the universe is 15 Gyr old (Jha et al. 1999). The final column 
of the table lists the contribution to the total black-hole merger rate from each cluster 
category. 



3.1. Merger rate in the local universe 

Given the black-hole merger rate corresponding to each category of star cluster, we 
now estimate the total merger rate TZ per unit volume. Table|2| lists, for various tyes of 
galaxies, the space densities and Sn, the specific number of globular clusters per M v = —15 
magnitude (van den Bergh 1995): 

S N = N GC 10 0A{M " +15) (4) 

(where Nqc is the total number of globular clusters in the galaxy under consideration). 
The values given for Sn in Table |^ are corrected for internal absorption; the absorbed 
component is estimated from observations in the far infrared. The estimated number 
density of globular clusters in the universe is 

GC = 8.4/i 3 Mpc" 3 . (5) 

A conservative estimate of the merger rate of black-hole binaries formed in globular clusters 
is obtained by assuming that globular clusters in other galaxies have characteristics similar 
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to those found in our own. The result is 

K GC = 5.4 x 10- 8 h 3 yr" 1 Mpc~ 3 . (6) 

Irregular galaxies, starburst galaxies, early type spirals and blue elliptical galaxies 
all contribute to the formation of young populous clusters. In the absence of firm 
measurements of the numbers of young populous clusters in other galaxies, we simply use 
the same values of Sn as for globular clusters. The space density of such clusters is then 
4>ypc — 3.5 h 3 Mpc -3 , and the black hole merger rate is 

TZypc = 2.1 x 10~ 8 /i 3 yr" 1 Mpc" 3 . (7) 

We find that galactic nuclei contribute negligibly to the total black hole merger rate. 

Based on the assumptions outlined above, our estimated total merger rate per unit 
volume of black-hole binaries is 

n = 7.5 x l(T 8 /i 3 yr" 1 Mpc~ 3 . (8) 

However, this may be a considerable underestimate of the true rate. First, as already 
mentioned, our assumed number (~ 10~ 4 iV) of ejected black-hole binaries is quite 
conservative. Second, the observed population of globular clusters naturally represents only 
those clusters that have survived until the present day. The study by Takahashi & Portegies 
Zwart (2000) indicates that ~ 50% of globular clusters dissolve in the tidal field of the 
parent galaxy within a few billion years of formation. We have therefore underestimated the 
total number of globular clusters, and hence the black-hole merger rate, by about a factor 
of two. Third, a very substantial underestimate stems from the assumption that the masses 
and radii of present-day globular clusters are representative of the initial population. When 
estimated initial parameters (Table [l], bottom row) are used, the total merger rate increases 
by a further factor of six. Taking all these effects into account, we obtain a net black-hole 
merger rate of 

K ~ 3 x 10- 7 h 3 yr" 1 Mpc~ 3 . (9) 

We note that this figure is significantly larger than the current best estimates of the 
neutron-star merger rate. 



3.2. LIGO observations 

The maximum distance within which LIGO-I can detect an inspiral event is estimated 
to be 



A^chirp 



5/6 



i? eff = 18Mpc I^PJ (10) 



(K. Thorne, private communication). Here, the "chirp" mass for a binary with component 
masses m 1 and m 2 is M cri i rp = (m 1 m 2 ) 3//5 /( m i + m 2 ) 1 ^. For neutron star inspiral, 
rrii — m 2 — 1-4 M Q , so M c hir P = 1.22 M Q , R e ^ = 21 Mpc, and we obtain the detection rate 
mentioned in the introduction. For black-hole binaries with m x = rri2 = m^ = 10 M Q , we 
find M c hirp = 8.71 M Q , R e g = 109 Mpc, and a LIGO-I detection rate of about 1.7 h 3 per 
year. For h ~ 0.65 (Jha 1999), this results in about one detection event every two years. 
LIGO-II should become operational by 2007, and is expected to have i? er j about ten times 
greater than LIGO-I, resulting in a detection rate 1000 times higher, or about one event 
per day. 



4. Discussion 

Black-hole binaries ejected from galactic nuclei, the most massive globular clusters 
(masses ^ 1O 6 M ), and globular clusters which experience core collapse soon after 
formation, tend to be very tightly bound, have high eccentricities and merge within a 
few million years of ejection. These mergers therefore trace the formation of dense stellar 
systems with a delay of a few Gyr (the typical time required to form and eject binaries), 
making these systems unlikely candidates for LIGO detections, as the majority merged long 
ago. This effect may reduce the current merger rate by an order of magnitude, although 
more sensitive future gravitational wave detectors may see some of these early universe 
events. We estimate that the most massive globular clusters contribute about 90% of the 
total black hole merger rate. However, while their black-hole binaries merge promptly upon 
ejection, the longer relaxation times of these clusters mean that binaries tend to be ejected 
much later than in lower mass systems. Consequently, we have retained these binaries in 
our final merger rate estimate (Eq. Q). But we note that represents a significant source if 
uncertainty. 

By the time the black hole binary is ejected it has experienced ~ 40-50 hard encounters 
with other black holes, as well as a similar number of encounters with other stars or 
binaries. During each of these latter encounters, there is a small probability that a low-mass 
star may collide with one of the black holes. Such collisions tend to soften the black hole 
binary somewhat (see Portegies Zwart et al. 1999), but they are unlikely to delay ejection 
significantly. A collision between a main-sequence star and a black hole may, however, lead 
to brief but intense X-ray phase. 

Finally, we have assumed that the mass of a stellar black hole is 10 M . Increasing this 
mass to 18M decreases the expected merger rate by about 50% — higher mass black holes 
tend to have wider orbits. However, the larger chirp mass increases the signal to noise, and 
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the distance to which such a merger can be observed increases by about 60% and the overall 
detection rate on Earth increases by about a factor of three. For 6 M black holes, the 
detection rate decreases by a similar factor. For black-hole binaries with component masses 
£ 12 M , the first generation of detectors will be more sensitive to the merger itself than to 
the inspiral phase that precedes it (Flanagan & Hughes 1998). Since the strongest signal is 
expected from black-hole binaries with high-mass components, it is critically important to 
improve our understanding of the merger waveform. Even for lower-mass black holes (with 
rribh £ 10 Mq), the inspiral signal comes from an epoch when the holes are so close together 
that the post-Newtonian expansions used to calculate the wave forms are unreliable (Brady 
et al. 1998). 
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Table 1: Overview of selected parameters for young populous clusters (Massey & Hunter 
1998), globular clusters (Djorgovski & Meylan, 1994) and galactic nuclei (McGinn et 
al. 1989). The first three columns list the cluster type, the total mass (in solar units) 
and the virial radius (in pc). The total mass and virial radius are given as distributions with 
a mean and the standard deviation around the mean. The orbital separation (in solar units) 
for a 1000 kT binary consisting of two 10 M black holes is given in the fourth column. The 
fifth and sixth columns list the expected number of black-hole binaries that are formed by 
the cluster and the fraction of these binaries which merge within 12 Gyr (allowing ~ 3 Gyr 
for the formation and ejection of the binaries and assuming a 15 Gyr old universe, Jha 1999). 
The contributions to the total black hole merger rate per star clusters per year (MR) are 
given in the final column. The bottom row contains estimated parameters for the zero-age 
population of globular clusters in the Galaxy, indicated with *. 
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~7 


< 


r^j O.O 


2500 


100% 


0.21 


Globular* 


6.0 ±0.5 


0±0.3 


33 


500 
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Table 2: Galaxy morphology class, space densities, average absolute magnitude (Heyl et 
al. 1997), and the specific frequency of globular clusters Sn (from van den Bergh, 1995 and 
McLaughlin 1999). The final column gives the contribution to the total number density of 
globular clusters. 
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